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Abstract: We present an updated assessment of the reach of the CERN LHC pp coUider 
for super symmetric matter in the context of the minimal supergravity (mSUGRA) model. 
In addition to previously examined channels, we also include signals with an isolated photon 
or with a leptonically decaying Z boson. For an integrated luminosity of 100 fb~^, values 



of m 



1/2 



1400 GeV can be probed for small mg, corresponding to a gluino mass of nig ~ 3 



TeV. For large mo, in the hyperbolic branch/focus point region, ~ 700 GeV can 

be probed, corresponding to rrig ~ 1800 GeV. We also map out parameter space regions 
preferred by the measured values of the dark matter relic density, the 6 ^ 57 decay rate, 
and the muon anomalous magnetic moment a^, and discuss how SUSY might reveal itself 
in these regions. We find the CERN LHC can probe the entire stau co-annihilation region 
and also most of the heavy Higgs annihilation funnel allowed by WMAP data, except for 
some range of large mo and m]^/2 if t^-ii ~ 50. 
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1. Introduction 



The search for supersymmetric (SUSY) matter is one of the primary objectives of experi- 
ments at high energy colliders Q. SUSY matter may reveal itself through indirect effects Q, 
as in contributions to rare decays such as 6 — > ,37, b — > sii or Bg — > , or via contribu- 

tions to electric or magnetic moments such as the dipole electric moment of the electron 
or neutron or (g — 2)^. It is possible that relic SUSY cold dark matter (CDM) has al- 
ready been detected gravitationally, and recent analyses of WMAP and other data sets 
indicate that the relic density f^cDAf^^ = 0-1126^Q'Qj;gj;(2(T CL)[^, with a baryonic den- 
sity about six times smaller. Both direct and indirect searches for relic SUSY CDM are 
underway Q. Detection of a signal in any of these experiments would provide evidence for 
physics beyond the standard model, but the results may have both supersymmetric as well 
as non-supersymmetric interpretations. 

The definitive discovery of SUSY matter will likely have to come from experiments 
operating at high energy colliders. Already, negative searches for SUSY by the LEP2 
experiments have resulted in significant bounds: for instance, the light chargino Wi must 
have mass > 103.5 GeV[P, while a SM-like higgs boson must have mass nih > 114.1 
GeV[^. The reach of the Fermilab Tevatron collider has been examined as well. In the 
context of the minimal supergravity model (mSUGRA)[^, with parameters 

mo, mi/2, ^0, tan/?, sign{n), (1.1) 

1711/2 values of up to 250 GeV (corresponding to rrig ~ 600 GeV) can be probed with 25 fb^^ 
for mo ~ 200 GeV(^ and small values of tan p. Then, the best reach is obtained via the 
clean trilepton channel. The reach is considerably reduced for larger values of mo or for 
large values of tan /?|9|, [l^, 11 1. 



The CERN LHC is expected to begin operating in 2008 with pp collisions at ^/s = 14 
TeV. While the initial luminosity is expected to be ~ 10 fb^^ per year, an integrated 
luminosity of several hundred fb~^ is ultimately anticipated. The reach of the LHC for 
supersymmetric matter has been evaluated in the mSUGRA model for low[12, 14, 15| 



and high[|T^, 13, 14| values of tan/3, for slepton[p!7|] and chargino-neutralino production ||l8|], 
and even in the experimentally unfavorable case where the lightest supersymmetric particle 
(LSP) decays hadronically via i?-parity violating interactions [p!9[ . The LHC reach has also 
been evaluated for models with non-universal soft terms as given by non-minimal SU (5) 
SUSY GUTs[20| and gaugino mediation models |pl[. Finally, the LHC reach has been 



evaluated for gauge-mediated [E^] (GMSB) and anomaly-mediated (AMSB) SUSY breaking 



models 12 3||. 



In studies evaluating the reach of the LHC, the signal channels have been classified by 
the number of isolated leptons present in each event. The isolated leptons usually arise as 



end products of cascade decays of gluinos, squarks or other massive SUSY particles 1 24 1. At 
large values of the parameter tan/3, the b and r Yukawa couplings become large, so that 
cascade decays to final states containing 6-jets and r-leptons are enhanced[25|. 



In this paper, we update our reach projections for the CERN LHC for several reasons. 
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We calculate the sparticle production and cascade decay events using ISAJET v7.64[26| 
This version includes a variety of radiative corrections and improvements to the 
sparticle mass spectrum that were not present in earlier ISAJET versions, including 
two-loop evolution of all RGEs. ISAJET 7.64 gives good agreement with spectra 
generated by the Suspect, SoftSUSY and Spheno codes, as compiled by Allanach, 
Kraml and Porod[^. In addition, 3-body decay matrix elements are included]^ 28|, 
so that decay product energy distributions are more accurately modeled. 

We adopt the code CMSJET v4.801 to model the CMS detector. This gives a more 
accurate portrayal of CMS than the toy detector models used in Ref. 

We include additional channels to our reach projections, including events containing 
a reconstructed Z — > l~ candidate = e or ^)p9[, and events containing isolated 



photons. The photonic events may arise from radiative neutralino decay pO|, 31 1 
Zi — > Zi7 which is enhanced in the low "hyperbolic branch" or "focus point" 

(HB/FP)[^, ^] region of parameter space. 

• Our parameter space scans extend over a wider range than earlier reach projections. 
We are motivated to do so to cover the HB/FP region with small [i where SUSY 
phenomenology can be significantly different. Further, we adopt a higher integrated 
luminosity value of 100 fb~^ than our earlier studies; this integrated luminosity should 
be achieved after several years of LHC operation. 

• We identify regions of mSUGRA parameter space consistent with recently updated 
constraints]^ from h — > 57, i^g — 'i)^ and '^^^^ calculations as well as the most recent 
constraints from the LEP2 experiments, and where there should be observable signals 
at the LHC. 

In Sec. 1^, we present the details of our computer calculations, cuts and detector 
simulation. In Sec. ^, we show our results for the LHC reach projections in the mo vs. mx/2 
plane, and identify parameter regions consistent with all exerimental constraints. We 
summarize our results in Sec. ^ 

2. Calculational details 

We use the CMSJET (version 4.801) fast MC package[||] for the CMS detector response 
simulation. Both SM background and signal events were generated using ISAJET 7.64 
which has undergone numerous upgrades since ISAJET 7.37 used in our last study [ 16 1. The 



improvements include the incorporation of matrix elements for the calculation of three body 
decays of gluinos, charginos and neutralinos, two loop renormalization group evolution of 
all couplings and soft SUSY breaking parameters, inclusion of 1-loop self energies for third 
generation fermions, and improved evaluation of ruA which, in turn, significantly moves 
the boundary of the allowed parameter space of the mSUGRA model. 

We have computed SM backgrounds from the following sources : tt, QCD 2 — > 2 
including cc and hb production, W + jets and Z + jets. Backgrounds from vector boson 
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pair production are negligible to the jetty signals that we consider in this study[12|. For 
the photonic signal we also included backgrounds from Wj and Z7 production, but found 
these to be negligible after hard cuts. Since the cross section of events with low pT is much 
larger than that of events with high pT, we generated the background events in several bins 
of pt- We followed the division described in |13]. We have also adopted the cuts suggested 



in this study for E^^^^, lepton, 1 lepton, 2 OS and SS lepton, 3 lepton and > 4 isolated 
lepton signals. We regard a lepton or a photon to be isolated if 

• it has no charged particle with pT > 2 GeV in a cone with AR < 0.3 around the 
direction of the lepton. 

• E£;ff" in the region with 0.05 < AR < 0.3 around the lepton's direction has to be 
less than 10% of the lepton transverse energy. 

For the convenience of the reader, we present these cuts in this paper. All events have 
to pass the following pre-cuts: 

• EJp^'' > 200 GeV; 

• Number of jets, Nj > 2. 

We use a modified UAl-jetfinder routine GLOBJF, implemented in CMSJET, to identify 
calorimeter jets. A cluster of particles is labeled as a jet if it has transverse momentum pT 
greater than 40 GeV and |r/| < 3. Leptons are required to satisfy the following pre-cuts: 

• Pt > W GeV for the muons, pT > 20 GeV for the electrons, \r]\ < 2.4 for both muons 
and electrons. 

• Electrons have to be isolated. Muon isolation is not required as part of the pre- 
cuts. We call an electron or muon non- isolated even if it satisfies the lepton isolation 
criteria but is part of a jet, or if it is isolated in the calorimeter, but non-isolated in 
the tracker. Naturally, if it is not isolated in the calorimeter then the lepton is called 
non-isolated. 

The pre-cuts for the photons are: 

• > 20 GeV in |r/| < 2.4. 

• photons have to be isolated. 

After events pass the pre-cuts, we impose 90% lepton detection efficiency for each lepton. 

The events which pass the pre-cuts are divided into signal types according to the 
number of leptons (or photons for the isolated 7 signal). In the case of the EJp'^^^ signal 
there can be any number of leptons, lepton signal has no leptons, 1 lepton signal has 1 
lepton, 2 OS lepton signal has 2 opposite sign leptons, 2 SS lepton signal has 2 same sign 
leptons, 3 lepton signal has 3 leptons, > 4 lepton signal has more than 3 leptons, Z — > 
signal has at least 2 OS, same flavor leptons with the invariant mass of this pair within the 
interval {Mz — AMz,Mz + AMz) [AMz is varied during the optimization procedure). 
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Finally, the isolated 7 signal has any number of leptons plus at least one photon (the cut on 
the number of photons is varied during the optimization procedure). Since muon isolation 
has not been included as part of the pre-cuts, if we impose the muon isolation during the 
optimization procedure, the number of events for some signal types can change. 

A signal in any channel is considered to be observable if after our optimization proce- 
dure described below, 

• the number of signal events 5 > 10 for an integrated luminosity of 100 fb~^, and 

• S > where B is the corresponding number of background events. 

We optimize the signal in each channel by imposing additional cuts. The set of cuts that 
we examined for this purpose are listed in Table |l[ For the optimization of the signal 
Z signal, we have an additional cut: the invariant mass of the pair of opposite sign 

same flavor lepton pair has to be in the interval (M^ — AMz,Mz + AMz) and AMz is 
taken to be 3, 6, 9, 30 GeV. For the case of the isolated photon signal, in addition to the 
optimization using the cuts in Table |l], we also vary the number of photons which can be 
1, 2, 3, 4, > 5. However, for large values of mo and mi/2, events with A^^'^° > 1 are generally 
too rare to pass the requirement of S" > 10. 

For each mSUGRA point that we analyze, we pass the various signals through each one 
of the complete set of cuts just discussed. If the signal satisfies our observability criterion 
for any one of these cut choices, we consider it to be observable.^ 

Table 1: The set of cuts that we have examined for the optimization of the SUSY signal. Except 
for the muon isolation, the numbers refer to the lower bound on the quantity listed in the first 
column. 



Variable(s) 


Values 




2, 3, 4, 10 


rjrp 


200, 300, 400, 1400 GeV 


E^2 


40, 150, 300, 400, 500, 600, 700, 800, 900, 1000 GeV 
40, 80, 200, 200, 300, 300, 400, 400, 500, 500 GeV 




0, 20 deg. 


Circularity 


0, 0.2 


fj, isolation 


on, off 



3. Results 

3.1 Reach of the LHC in various channels 

Since sparticle masses are largely determined by the parameters mo and m^/2 the mo — 

^When more than one cut choices lead to an observable signal, we retain the choice that maximizes the 
quantity S/VS + B. 
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plane is a convenient arena for simultaneously displaying the SUSY reach of future exper- 
imental facilities together with regions already excluded by current data. Many previous 
calculations for the LHC reach are, therefore, presented in this plane (for various choices of 
other parameters) starting with tan /? = 2; this low of a tan (3 value is now largely excluded 
by the LEP2 bound on m/^. We begin our presentation with tan /? = 10, in Fig. |l|. The red- 
shaded region on the left is excluded because either electroweak symmetry is not properly 
broken, or fi is the LSP, while that on the right is excluded because radiative electroweak 
symmetry breaking (REWSB) does not occur. The magenta region in the lower left is 
excluded by LEP2 bounds on > 103.5 GeV, and nih > 114.1 GeV (for a SM-hke hght 
Higgs h). The maximum reach is shown by the E^^^^ contour, where the signal events in- 
clude all isolated lepton possibilities. Also shown for reference are contours of nig = 2 TeV, 
and rriq = 2 TeV. The maximum reach in mi/2 occurs for low mo, where squark masses are 
somewhat lighter than gluinos, so that qq, qg and gg production processes all have large 
rates. Gluinos as heavy as ~ 3 TeV may be detectable at the LHC if squarks are somewhat 
lighter. As itlq increases, squark masses increase, so gluino pair production becomes the 
dominant sparticle production mechanism. At very large values of mo, the squarks (and 
also sleptons) essentially decouple, and the reach contours flatten out, since rrig is roughly 
constant for each value of mi/2- The maximal reach for 100 fb~^ of integrated luminosity 
is mi/2 ~ 700 GeV, corresponding to a gluino mass value of rrig ~ 1800 GeV. 

The reach for SUSY signals in the individual channels introduced in the last section 
are shown by the various contours labelled by the corresponding topology in Fig. [l|. For 
low values of mo, sleptons and squarks are relatively light, and gg, gq and qq production 
processes all occur at large rates; direct slepton and sneutrino pair production rates are 
much smaller. However, left-squarks qi decay frequently into the heavier chargino and 
heavier neutralinos, which in turn may decay to sleptons. These cascade decays frequently 
terminate in isolated leptons, which together with leptons from decays from tops and stops 
in SUSY processes, result in large rates for leptonic signals. For a fixed value of mi/2, as mo 
increases, g decay to Ijl becomes suppressed, and BF{g qn) increases. Since qn — > qZi 
most of the time, cascade decays then give rise to a higher fraction of + jets events, 
and the leptonic signal from the decay of gluinos is reduced. Furthermore, since squarks 
become heavier relative to gluinos as mo is increased, the leptonic signal from directly 
produced qi also becomes smaller. As mo increases even further, g — > q^q also becomes 
suppressed or even forbidden, and g — > tit (and possibly g 661,2) dominates. The decays 
through tops and stops gives rise again to leptonic states due to ti and t leptonic decays. 
This results in an increased reach at moderate mo values via leptonic modes such as SS 
and OS dileptons, and 3£ events. As mo increases even more, g two-body decays become 
completely forbidden, and three-body decays dominate, and the leptonic reach contours 
tend to level off. 

Finally, we show a contour that marks the signal reach for events including + 
jets plus at least one isolated photon. In evaluating this reach, we have only retained 
physics backgrounds in our calculation. Detector-dependent backgrounds where a jet fakes 
a photon may be significant.^ In most of the parameter space, the additional photon arises 
^If we take the probability for a jet to fake a photon to be 5 x 10"'* and assume that the hard scattering 
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mSugra with tan^ = 10, = 0, [i> 




rUf. ( GeV) 



Figure 1: The reach of CERN LHC in the mo vs. toi/2 parameter plane of the mSUGRA model, 
with tan/3 — 10, Aq — and /i > 0, assuming 100 fb~^ of integrated luminosity. The red (magenta) 
regions are excluded by theoretical (experimental) constraints discussed in the text. We show the 
reach in the 0£, 1£, OS, SS, 31, > 4£, 7 and Z channels, as well as in the "inclusive" channel. 

from /i ^ 77 decay, where the h is produced copiously in sparticle cascade decays, especially 
from Z2 — > Zih. In these regions, in fact, if we require two isolated photons, then we can 
reconstruct a di-photon invariant mass. This is illustrated in Fig. ^ for the parameter space 
point ruQ = = 500 GeV, Aq = 0, tan/3 = 30 and /i > 0. It is amusing to note that the 
/i — > 77 signal should be visible in SUSY events for 100 fb~^ of integrated luminosity. We 
see that (for these parameters) the highest possible luminosity is needed for the detection 

events have ~ 10 30-40 GeV "jets" in them, about one in 500 background events will also appear to have 
an isolated photon. Assuming that this fake photon background can be estimated by reducing the physics 
background in the inclusive E™"" channel by 500, we find that this background is somewhat smaller, but 
of the same order of magnitude as the physics background that we have evaluated. A real evaluation of 
this detector-dependent background is beyond the scope of our analysis. 
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mSugra with mg=500 GeV, niy^^SOO GeV, tan^ = 30, = 0, |i > 




60 70 80 90 100 110 120 130 140 150 

m(YY) (GeV) 



Figure 2: The diphoton invariant mass in SUSY events witli two isolated photons. The shaded 
histogram is for diphotons from SUSY events, while the hatched histogram represents corresponding 
events of SM origin. The background has been obtained by scaling a Monte Carlo run for a lower 
value of integrated lumnosity. 



of this severely rate-limited, but essentially SM background-free, signal. A detection of h 
in this manner would nicely confirm its detection via gg ^ h ^ 77, where the signal would 
be picked out as a small peak above an enormous continuum background. A Higgs signal 
in the multijet plus 77 + E^^^^ channel would also suggest the supersymmetric origin 
of the Higgs production process. As niQ becomes very large, and the region of low is 
approached, the radiative decay Z2 — > Zi^ becomes enhanced|31]. In this HB/FP region, 
the isolated photon contour likewise turns up, to reflect the increase in isolated photon 
activity from neutralino radiative decay. 

As we move to larger values of the parameter tan /3, as in Fig. ^ for tan /3 = 30, the first 
thing to notice is that the left red-shaded region, where fi is the LSP, has expanded. This is 
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because as tan/? increases, both the r (and also b) Yukawa couphngs become non-neghgible. 
This results in a reduction of stau and sbottom soft breaking masses via RGE running, 
and also in greater L — R mixing, which again reduces the mass of the lightest eigenstates. 
The negative results of searches for heavy isotopes of hydrogen (or other elements) results 
in limits that are many orders of magnitude below their expected relic density in big-bang 
cosmology, so that this possibility is strongly excluded. Finally, the large values of fh and 
fr together with the reduction of the corresponding sfermion masses relative to the first 
two generations mentioned above, increases various three body sparticle decays to 6-quarks 
and T leptons, at the expense of their first and second generation counterparts ||2^. The 
net effect of this is to enhance sparticle cascade decays to final states containing 6-quarks 
and T-leptons in the low mo region. This results in a diminution of the reach in isolated 
lepton channels at low mo, as compared with Fig. ||. The large h and r Yukawa couplings 
hardly affect the + jets and 0£ signals, and so the ultimate reach of the LHC changes 
little in proceeding from tan/3 = 10 to tan/3 = 30. 

In Fig. ^, we show the mo vs. ^1^12 plane for tan/? = 45 and /x < 0. The large h 
and r Yukawa coupling effects are accentuated even more in this figure: a larger region 
is excluded at low mo, and the reach via multi-lepton channels is further diminished for 
small mo values where squarks and sleptons still play a role in determining cascade decay 
patterns. However, again, the overall reach in the + jets and 0£ channels is hardly 
affected. If we increase tan /? much beyond about 50 for /i < 0, the parameter space begins 
to close up fast due to a breakdown in REWSB. 

In Fig. m, we show the mo vs. mj^/2 plane for tan/3 = 52 and fi > 0. In this case, 
the overall reach in the ^Jt + jets and Oi channels is similar to the cases at lower tan (3. 
However, in the multi-lepton channels, there is again a suppression of reach at low mo. This 
is because for very high tan /3, m^j is so light that Wi dominantly decays to fiVr rather 
than ZiW, and likewise, Z2 — > fir rather that Zih or ZiZ. As mo increases, the stau mass 
increases, and the Wi — > fiUr decay mode becomes more suppressed, which increases the 
Wi ZiW branching fraction. The subsequent W boson decays from Wi Z\W lead 
to hard isolated leptons. 

Our projections of the LHC reach, where they can be directly compared, are qualita- 
tively similar to the results in Ref. Differences between the results can be attributed 
to the difference in the observability criteria (S > h\J S + B used in Ref.[^] requires a 
minimum of 25 events to be compared with 10 events in our study, as well as a somewhat 
larger significance of the signal) , and to the use of PYTHIA instead of IS A JET for event 
simulation. 

3.2 The LHC reach in light of indirect constraints 

A variety of low energy measurements have been used to obtain constraints on the param- 
eter space of the mSUGRA model. These include the measured values of the cold dark 
matter density, the branching fraction BF(b S7), the value of the anomalous magnetic 
moment of the muon, = {g — 2)^/2, ^ and a lower limit on BF{Bs n~^fi~). Unlike 

^There is considerable theoretical uncertainty in the SM value of a^, so that caution must be exercised 
in interpreting the result of experiment E821. 
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Figure 3: The reach of CERN LHC in the mo vs. toi/2 parameter plane of the mSUGRA model, 
with tan/3 — 30, Aq = and /i > 0, assuming 100 fb~^ of integrated luminosity. The red (magenta) 
regions are excluded by theoretical (experimental) constraints discussed in the text. We show the 
reach in the 0£, 1£, OS, SS, 31, > 4£, 7 and Z channels, as well as in the "inclusive" channel. 



collider constraints which are much more direct, constraints from low energy measurements 
may be considerably more sensitive to details of the model, or to improvements in the the- 
oretical calculation. Within a specific framework {e.g. mSUGRA), however, these indirect 
constraints exclude certain regions of parameter space, and also suggest other regions where 
future searches might be focussed. 
Neutralino relic density: 

Measurements of galactic rotation curves, binding of galactic clusters, and the large 
scale structure of the universe all point to the need for significant amounts of cold dark 
matter (CDM) in the universe. In addition, recent measurements of the power spectrum 
of the cosmic microwave background from WMAP and other data setsQ lead to 

. ncDMh^ = o.n26t?,-Z9- 
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Figure 4: The reach of CERN LHC in the mo vs. toi/2 parameter plane of the mSUGRA model, 
with tan/3 = 45, Aq = and /i < 0, assuming 100 fb~^ of integrated luminosity. The red (magenta) 
regions are excluded by theoretical (experimental) constraints discussed in the text. We show the 
reach in the 0£, 1£, OS, SS, 31, > 4£, 7 and Z channels, as well as in the "inclusive" channel. 



The upper limit derived from this is a true constraint on any stable relic from the Big Bang, 
such as the lightest neutralino of the mSUGRA model. Regions of mSUGRA parameter 
space which result in a relic density that violates this bound are excluded. A remarkable 
feature of the model is that there are regions of the parameter space where the neutralino 
density lies in the observed range, so that the neutralino makes up almost all the cold 
dark matter in the universe. We remark, however, that unlike the upper limit above, the 
corresponding lower limit is flexible, since there may be additional sources of CDM such 
as axions, or states associated with the hidden sector of the mSUGRA model and/or extra 
dimensions. 

To estimate the relic density of neutralinos in the mSUGRA model, we use the recent 
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Figure 5: The reach of CERN LHC in the mo vs. toi/2 parameter plane of the mSUGRA model, 
with tan/3 — 52, = and /i > 0, assuming 100 fb~^ of integrated luminosity. The red (magenta) 
regions are excluded by theoretical (experimental) constraints discussed in the text. We show the 
reach in the 0£, 1£, OS, SS, 31, > 4£, 7 and Z channels, as well as in the "inclusive" channel. 



calculation in Ref. [35|. All relevant neutralino annihilation and co-annihilation reactions 
are included along with relativistic thermal averaging 1 36], which is important for obtain- 
ing the correct neutralino relic density in the vicinity of annihilations through s-channel 
resonances. 

BF(6^ S7): 



The branching fraction BF[h sj) has recently been measured by the BELLE|37|, 
CLEO|^^ and ALEPH||3^ collaborations. Combining statistical and systematic errors in 
quadrature, these measurements give (3.36 ± 0.67) x 10"'' (BELLE), (3.21 ± 0.51) x lO"'' 
(CLEO) and (3.11 ± 1.07) x 10"^ (ALEPH). A weighted averaging of these results yields 
BF{b 57) = (3.25 lb 0.37) x 10~^. The 95% CL range corresponds to ±2(7 away from 
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the mean. To this we should add uncertainty in the theoretical evaluation, which within 
the SM dominantly comes from the scale uncertainty, and is about 10%. Together, these 
imply the bounds, 

• 2.16 X 10"^ < BF{b S7) < 4.34 x 10-"^. 

The evaluation of the SUSY contribution to this decay entails additional theoretical un- 
certainties, especially when tan/3 is large, so that this range should be relaxed somewhat. 
In our study, we show contours of BF{b — > 57) of 2, 3, 4 and 5 x 10~^. 

The calculation of BF{h S7) used here is based upon the program of Ref. In 
our calculations, we also implement the running 6-quark mass including SUSY threshold 
corrections as calculated in ISAJET; these effects can be important at large values of the 
parameter tan/3[^]. Our value of the SM 6 — > 57 branching fraction yields 3.4 x lO"'^, 
with a scale uncertainty of 10%. 

Muon anomalous magnetic moment 

The muon anomalous magnetic moment = [g — 2)^/2 has been recently measured 
to high precision by the E821 experiment Q: = 11659204(7) (5) x 10"^°. The most 
challenging parts of the SM calculation are the hadronic light-by-light 1 43] and vacuum 



polarization (IIVP)[^] contributions and their uncertainties. Presently these results are 
in dispute. In the case of the HVP the use of tau decay data can reduce the error, but 
the interpretation of these data is somewhat controversial^]. Thus, the deviation of the 
measurement from the SM depends on which prediction is taken into account. According 
to the recent analysis by Hagiwara et al. [^ : 

• 11.5 < 5a^ X 10^0 < 60.7. 

A different assessment of the theoretical uncertainties! 44] using the procedure described in 
Ref.|4|] gives, 

• -16.7 < 5a^ X 10^° < 49.1. 

Yet another determination has recently been made by Narison, includes additional scalar 



meson loopsl47]. This, using e"^e — > hadrons data to evaluate hadronic vacuum polariza- 



tion contributions, yields 

• -3.9 < 6a^, X 10^° < 52.1, 
while using r-decay data results in 

• -26.7 < 5a^ X 10^° < 30.1. 

The latter may include additional systematic uncertainties from how isospin breaking effects 
are incorporated. 

In view of the theoretical uncertainty, we only present contours of (5a^, as calculated 
using the program developed in [^], and leave it to the reader to decide the extent of the 
parameter region allowed by the data. 

Eg — > /U^^^ decay 

The branching fraction of Bg to a pair of muons has been experimentally bounded by 
CDF[||]: 

• BF{Bs ^+/u") < 2.6 X 10"^ 
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If tan /3 ~ 20 — 25 SUSY contributions to this decay are small and do not lead to new 
constraints on the parameter space. If tan f3 is large, the important SUSY contribution 
to this decay is mediated by the neutral states in the Higgs sector of supersymmetric 
models. While this branching fraction is very small within the SM {BFsm{Bs /^^/^ ) — 
3.4 X 10~^), the amplitude for the Higgs-mediated decay of roughly grows as tan^ /?, ^ 
and hence can completely dominate the SM contribution if tan /? is large. In our analysis 
we use the results from the last paper in Ref. [^] to delineate the region of mSUGRA 
parameters excluded by the CDF upper limit on its branching fraction. 

In Fig. ^ we again show the mo vs. mi/2 plot for tan/3 = 10, Aq = and ^ > 0. 
This time, we exhibit contours for the low energy observables mentioned above, as mapped 
out in Ref. [^^. A analysis of the indirect constraints was performed in Ref. pT| , 
which helped to identify regions of parameter space allowed by all the combined indirect 
constraints. Several regions emerge in the niQ vs. mi/2 plane where the relic density can 
be within the WMAP range. At low mo and low mi/2, neutralino annihilation through 
t-channel sleptons can occur with a high rate. Much of this so-called "bulk" region is 
largely excluded by the LEP2 Higgs bound (shown as the red contour), and in addition, 
BF{h 57), though in the acceptable range, is below its experimental central value. In 
this region, sparticles are very light, and a SUSY discovery by the CERN LHC should be 
easy. 

One of the remaining regions allowed by relic density constraint is the very narrow strip 
adjacent to the "Zi not LSP" region at low mo, where stau co-annihilation is important |^2|| . 
We see that the reach of the LHC apparently covers much of the stau co-annihilation strip, 
up to mi/2 ~ 1400 GeV for 100 fb^^ of integrated luminosity. To determine the LHC reach 
needed to completely explore the stau co-annihilation strip, we show the relic density for 
fixed mi/2 values vs. mo in Fig. |^. In frame a), it is evident that the stau co-annihilation 
corridor yields a relic density in accord with WMAP results for mi/2 ~ 900 GeV. Thus, 
the LHC would completely explore the stau co-annihilation corridor for tan f3 = 10. 

The other region allowed by relic density constraint occurs at very large mo where 
/i becomes small: the HB/FP region |33t|.^ Here, the growing higgsino component of 



the neutralino allows for efficient annihilation into vector boson pairs 1 53, For very 



small fi values, then neutralino-chargino co-annihilation becomes important [pqj. In this 
FP/HB region squarks are very heavy, and as mi/2 increases, nig increases as well so 
that strongly interacting sparticle production cross sections decrease. Since fi is small, the 
light charginos and neutralinos are significantly higgsino-like, and can only be produced 
via electroweak interactions. We see that the reach of the LHC is limited to mi/2 ~ 700 
GeV. For even higher mi/2 values, and staying in the HB/FP region, one enters an area 
that is not accessible to LHC experiments, at least via the general purpose search strategies 
described here, even though the light charginos and neutralinos are expected to have masses 
less than 250-500 GeV. It would be interesting to examine whether it is possible to devise 



''The tan"^ {3 growth obtains if the tree-level value of mt is fixed. For large values of tan l3, the radiative 
correction to rrib is important resulting in a deviation from the tan"^ P growth. 

^The entire HB/FP region is not shown here: over part of this region, the evaluation of the fi parameter 
using ISAJET 7.64 is numerically unstable. This is corrected in ISAJET v7.65. 



- 13 - 



mSugra with tan^ = 10, A^ = 0, \i> 




1000 2000 3000 4000 5000 



Mg ( GeV) 

— m^=114.1GeV — a^susYxlO^^ — Br(b^sy)xlO^ 
Q.h^= 0.094 0.129 — 1.0 

Figure 6: Contours of several low energy observables in the toq vs. mi/2 plane of the mSUGRA 
model, for tan (3 — 10, Aq — and /i > 0. We show contours of CDM relic density (green), together 
with a contour of irih = 114.1 GeV (red), contours of muon anomalous magnetic moment (x lO^'') 
(blue) and contours of 6 ^ 57 branching fraction (xlO'*) (magenta). Also shown is the maximal 
reach of the CERN LHC in the ^r+jets channel for 100 fb^^ of integrated luminosity. 

search strategies that exploit specific characteristics of the FP/HB region to extend the 
LHC reach. For instance, in the small region, decays of gluinos into third generation 
fermions are enhanced resulting in events with hard b and t quarks from the decay of heavy 
gluinos. Whether it is possible to extend the LHC reach in these event topologies merits 
further investigation. 

Next, we show the contours for the same low energy observables as in Fig. ^ but for 
tan /3 = 30, ^0 = and > in Fig. |^. The larger value of tan /? causes the values 
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Figure 7: The relic density ^^_^h^ vs. mg for various fixed mi/2 values in the stau co-annihilation 
corridor, for a) tan/3 — 10, jj. > 0, b) tan/3 — 30, /i > 0, c) tan/3 — 45, /i < and tan/3 — 52, /i > 0. 
We also show via green lines the 2a WMAP limits on il^^ /i^ . 



and BF[b 57) to deviate more from their SM values at low niQ and mi/2', the bulk 
region of relic density annihilation is disfavored by BF{h — > S7). This leaves only the stau 
co-annihilation corridor and the HB/FP region as phenomenologically viable. The LHC 
reach contour remains similar to that for tan/3 = 10. Again, we see that if SUSY has 
parameters in this plane, then LHC should see it unless > '''00 GeV with SUSY in 
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Figure 8: The same as Fig. o but for tan/3 = 30. 



the HB/FP region. In particular, the LHC should be able to cover the stau co-annihilation 
corridor for tan/3 = 30, since it yields a relic density ^^J^'^ ^ 0.129 for ~ 1050 GeV, 
while the LHC reach for low uiq extends to mi/2 ~ 1400 GeV. 

Fig. P illustrates the impact of the low energy observables, this time for tan/3 = 45 
and ^ < 0. The LHC reach is similar to the lower tan [3 cases, but now the low mo and 
m-x/2 bulk region is firmly excluded by both and BF(b sj), which are large negative, 
and large positive, respectively. We see, however, that a new region of low relic density 
has opened up: the ZiZi ^ A, H ^ ff annihilation corridor, where annihilation takes 



place especially through the very broad A width|53, 54]. In fact, we see that although the 
A annihilation corridor extends to very large values, the region allowed by WMAP 
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Figure 9: The same as Fig. ^ but for tan /3 = 45 and /i < 0. 



is almost entirely accessible to LHC searches. A modest additional integrated luminosity 
beyond 100 fb~^ assumed in this study should cover this entire region. In addition, the stau 
co-annihilation corridor remains consistent with WMAP constraints up to mi/2 values as 
high as 1200 GeV as shown in frame c) of Fig. 0, so that LHC should be able to completely 
explore this region. Finally, the HB/FP region can again be explored up to mi/2 ~ 700 
GeV at the LHC. 

we show the same low energy contours in the mo vs. mi/2 plane, but 
52 and ^ > 0. In this case, the low mo and mx/2 bulk region is largely 
excluded because BF{b — > S7) < 2 x 10~^. The WMAP constraints again restrict us to 
either the stau co-annihilation region at low mo, or the HB/FP region at large mo. We 



In Fig 
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can see from Fig. that the stau co-annihilation corridor extends to mi/2 ~ 1870 GeV, 
which is somewhat beyond the reach of the LHC for 100 fb~^ of integrated luminosity. 
Models with these high parameter values suffer from considerable fine-tuning, since the ^ 
parameter, which has been invoked as a measure of fine-tuning [^2|, is beyond 1600 GeV 
for mi/2 > 1400 GeV, so that /M"^ is large. In the HB/FP region, the LHC reach, via 
these channels, again cuts off at mi/2 ~ 700 GeV: once again, it would be worth exploring 
whether the SUSY signal can be picked up in other channels. If we increase tan (3 beyond 



52, then the A, i?-annihilation funnel re-enters the figure [pl| . In Ref. |51|, it is shown 
that this annihilation funnel can reach values as high as 1400 GeV for mo ~ 3.5 TeV 
and tan/3 = 56. Thus, in this case, the LHC will not be able to access the complete A, H 
annihilation funnel. 



4. Conclusions 

We have updated our assessment of the SUSY reach of the CERN LHC via EJp'^'^^ and 
multilepton channels, and presented new results for the reach in channels with isolated 
photons or leptonically decaying Z bosons. We work within the framework of the mSUGRA 
model, and use ISAJET v7.64 together with the CMSJET fast detector simulation to 
model the CMS detector, and assume an integrated luminosity of 100 fb~^. Our results are 
presented over an expanded mSUGRA model parameter space to include the reach in the 
so-called HB/FP region at very large mg. This region, together with the stau coannihilation 
corridor, and the annihilation funnel where LSPs annihilate via the Aoi H resonances, are 
strongly preferred by the recent data of the WMAP collaboration. The overall LHC reach 
turns out to be quite insensitive to tan (3. We find that experiments at the LHC will probe 
mi/2 ~ 1400 GeV for low mg, and mij2 ~ 700 for large mo in the HB/FP region. These 
values correspond to rrig ~ 3000 GeV and 1800 GeV, respectively. 

We have also presented the reach in a variety of multi-lepton channels. The reach in 
these individual channels is, in general, sensitive to tan (3. We also show the reach in a 
channel including reconstructed U decays, and channels including isolated photons. 

The isolated photon signals may contain /i — > 77 events at a low, but observable rate. 
Indeed the SUSY event sample may contain SM background-free /i ^ 77 events, though 
a very high integrated luminosity will be needed to identify the rate-limited h signal. In 
the HB/FP region, the photonic channels also have a slight enhancement from radiative 
neutralino decays Z2 — > ^17. For mi/2 ~ 800(400) GeV and small (large) values of mo, 
there should be observable signals in all these channels if new physics discovered at the 
LHC is to be interpreted as supersymmetry as realized in the mSUGRA model. 

We have also examined the reach of the LHC in light of the recent assessment of direct 
(from LEP2) and indirect constraints on the mSUGRA model. The indirect constraints 
include the neutralino relic density ^^^h'^ from recent WMAP analyses together with 
accelerator measurements of BF{h — > 57), = {g—2)^ and the bound on BF{Bs fJ-^fJ- ) 
(this bound is hardly constraining for the parameter planes that we have examined). For 
large values of tan /3, experimental values of BF{b S7) and disfavor negative values 
of the fi parameter unless mo and mi/2 are also large, but for tan/3 ~ 10, values of mo 
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Figure 10: The same as Fig. g but for tan/3 = 45. 



and mi/2 ~ 400 — 500 GeV are perfectly acceptable. For /x > 0, mo would have to be 
rather small so that the relic density is either in the bulk annihilation region or in the stau 
coannihilation strip, or niQ would have to be very large, in the HB/FP region. The CERN 
LHC can definitely explore all the bulk annihilation region, and can explore all the stau 
co-annihilation corridor unless tan /3 is very high. The LHC can explore the HB/FP region 
up to mi/2 ~ 700 GeV via the conventional SUSY search channels. However, the HB/FP 
region appears to extend indefinitely to large mi/2 and niQ values, ultimately well beyond 
the LHC reach. For large tan the A, H annihilation funnel enters the mo vs. plane. 
The LHC with 100 fb"^ ought to be able to detect SUSY over this entire region, except 
for the case of very large tan/3 ~ 56 with /i > 0. 
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